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Digital  10  cm  radar  data  were  collected  using  the  Texas  A&M 
University  weather  radar  system.  The  digital  radar  data  were  then 
compared  with  GOl'S-East  infrared  imagery  mapped  to  the  digital  radar 
grid  projection.  Point  values  of  infrared  digital  count  were  correlated 
with  point  values  of  zero-tilt  reflectivity,  low-level  PVSZ,  mid-level 
PVSZ,  upper-level  PVSZ,  and  VIL  for  selected  portions  of  two  digital 
radar  tilt-sequences.  Correlations  also  were  performed  using 
thresholded  infrared  images  and  thrcsholded  digital  radar  data. 

The  results  indicated  a  positive  correlation  between  infrared 
digital  count  and  the  digital  radar  data.  Thresholding  the  infrared 
Imagery  was  found  to  improve  the  correlation  between  infrared  digital 
count  and  digital  radar  data.  Thresholding  the  digital  radar  data  pro¬ 
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CHAPTER  I 


INTRODUCTION 

The  importance  and  use  of  weather  radar  as  a  means  of  detect¬ 
ing,  observing,  and  forecasting  severe  weather  phenomena  is 
firmly  established.  Meteorological  satellite  data ,  particularly 
data  from  geostationary  satellites,  have  also  proved  to  be  of  groat 
value  in  detecting,  observing,  and  forecasting  severe  weather 
phono  lena.  The  increasing  availability  of  digital  data  from  both 
r.rjtc oroiogical  radars  and  geostationary  satellites  makes  ready  com- 
pan  on  of  both  types  of  data  feasible  and  useful  in  the  analysis  of 
severe  weather  phenomena.  The  possibility  that  characteristics  of 
severe  weather  phenomena  shared  in  common  by  the  two  types  of 
data  can  be  used  in  the  real-time  analysis  and  forecasting  of  severe 
weather  has  made  comparisons  of  these  data  attractive  to  both 
operational  and  research  meteorologists. 

1.  The  Need  for  This  Investigation 

The  increasing  automation  of  systems  for  handling  large  volumes 
of  meteorological  data  and  the  proliferation  of  interactive  systems 

The  citations  on  the  following  pages  follow  the  style  of  the 
Journal  of  Applied  Meteorology. 
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designed  to  display  and  compare  fields  of  satellite,  digital  radar, 
and  conventional  surface  and  upper-air  data  have  focused  increased 
interest  on  the  development  of  procedures  to  acquire  and  analyze 
such  data  on  a  real-time  basis  for  the  forecasting  of  severe  weather 
development.  Interactive  meteorological  data  handling  systems  are 
in  use  at  some  universities  and  within  government  agencies  charged 
with  research  and  development  responsibilities.  Interactive  systems 
are  expected  to  be  in  use  at  selected  central  forecasting  facilities 
within  the  National  Weather  Service  and  the  Air  Weather  Service  in 
the  near  future.  These  systems  must  allow  for  rapid  analysis  of  all 
data  fields  and  ready  comparison  of  significant  features  in  each 
field  of  data  with  features  in  the  other  data  fields.  Maximizing  the 
usefulness  of  such  systems  in  forecasting  applications  also  requires 
an  understanding  of  the  significant  relationships  between  data  fields 
such  as  digital  satellite  data  and  digital  radar  data.  The  large 
volume  of  high  resolution  digital  data  available  at  frequent  intervals 
from  the  geostationary  satellites  such  as  the  GOER  (Geostationary 
Operational  Environmental  Satellite)  spacecraft  suggests  that  signi¬ 
ficant  relationships  may  be  established  between  the  satellite  data 
and  simultaneous  digital  radar  data  which  may  be  useful  in  the 


development  of  rainfall  and  flood  forecasting  techniques. 


2.  Present  Status  of  Research  Relating  to  This  Investigation 

Comparisons  of  meteorological  satellite  imagery  with  radar  data 
have  generally  been  either  of  a  qualitative  or  a  quantitative  nature. 
Qualitative  comparisons  began  shortly  after  the  launch  of  TIROS  1 
(Television  and  Infrared  Observation  Satellite  1),  the  world's  first 
full-time  meteorological  satellite.  Blockmer  (1961)  and  Nagle  (1903) 
used  TIROS  1  imagery  in  comparison  with  available  radar  precipita¬ 
tion  patterns  to  study  squall  line  activity.  Nagle  and  Serebreny 
(1962)  compared  TIROS  1  imagery  with  scope  photographs  from  ground- 
based  and  airborne  radars  in  studying  the  precipitation  patterns 
associated  with  a  maritime  cyclone.  Boucher  (1967)  compared  ima¬ 
gery  from  several  TIROS  spacecraft  with  available  radar  data  and 
mesoscale  analyses  in  examining  the  relationship  of  satellite-ob¬ 
served  cirrus  shields  to  the  severity  of  thunderstorm  systems. 

Boucher  found  that  the  diameter  of  the  cirrus  shield  resulting  from 
the  integration  of  all  anvil  clouds  within  a  thunderstorm  complex  was 
an  indicator  of  storm  severity.  Blackmer  and  Serebreny  (1966) 
developed  models  of  cloud  and  precipitation  patterns  associated  with 
maritime  cyclones  based  on  TIROS  9  imagery,  shipborne  and  shore- 
based  radar  data,  and  surface  and  upper-air  data. 

Vonder  Haar  (I960)  employed  data  from  the  geostationary  satel¬ 
lites  ATS  1  (Applications  Technology  Satellite  1)  and  ATS  3  in  a  study 


of  reflected  re  die  nee  explications  which  presented  a  detailed  com¬ 
posite  of  a  reflected  radiance  field  and  meteorological  radai  reports. 
Vondcr  Haar  and  Cram  (1970)  also  used  ATS  3  imagery  in  a  short- 
range  forecasting  study  which  included  radar  data  as  well  as  surface 
observations  within  a  mesonetwork.  SMS/GOL’S  (Synchronous 
Meteorological  Satollitc/GOES)  imagery,  rningage  data,  and  radar 
data  provided  the  data  base  for  Fujita's  (1977)  study  of  meteorologi¬ 
cal  satellite  applications  for  army  operations. 

Quantitative  comparisons  of  meteorological  radar  and  meteoro¬ 
logical  satellite  data  have  been  made  using  data  from  both  polar- 
orbiting  and  geostationary  satellites.  Woxlor  and  Allisor.  (1  972) 
employed  Nimbus  4  infrared  data  and  a  reflectivity  map  from  the 
National  Severe  Storms  Laboratory  (NSSL)  WSR-57  radar  to  study  a 
tornado  outbreak  that  occurred  near  Oklahoma  City.  Gruber  (1973) 
also  used  Nimbus  4  infrared  data  to  determine  areas  of  deep  convec¬ 
tion  and  compared  his  findings  to  rariem  data  as  part  of  a  study  con¬ 
cerned  with  estimating  rainfall  in  convective  areas. 

ESSA  5  (Environmental  Science  Services  Administration  5)  and 
ESSA  7  satellite  digital  cloud  mosaic  films  wore  examined  by  Gerrish 
(1970)  to  determine  the  cloudiness  in  each  element  of  a  grid  array 
centered  south  of  Miami,  Florida.  Gerrish  used  those  data  along 
with  WSR-57  radar  data,  raingage  data,  and  upper-air  data  in 
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examining  the  relationship  between  rainfall  and  cloudiness  within 
the  grid  area.  Gerrish  (1975)  later  used  DMSP  (Defense  Meteoro¬ 
logical  Satellite  Program)  satellite  visible  and  infrared  data  and 
ATS  3  visible  data  in  a  similar  study  of  cloudiness  and  precipitation. 
Gerrish  determined  that  radar  echo  area  correlated  best  with  DMSP 
infrared  cloud  area  at  18,000  ft  and  at  30,000  ft  and  that,  on  the 
average,  radar  echoes  occupied  only  42  per  rent  of  the  cloud  area 
at  those  altitudes.  DMSP  data  also  were  used  by  Phillips  (1  975) 
in  a  study  that  compared  infrared  satellite  data  with  digital  radar 
data  from  the  NSSL.  Phillips  concluded  that  the  DMSP  imagery  could 
be  studied  in  greater  detai  1  when  digitised  using  a  :v,  i  crod  en  silo  mot  r 
but  that  detailed  determination  of  the  structure  of  convective  storm 
activity  did  not  seem  possible  with  the  infrared  sensors  then  in  use. 

Blackmcr  (1975)  studied  the  extent  to  which  a  combination  of 
visible  and  infrared  data  indicated  the  presence  and  intensity  of 
precipitation  within  a  cloud  and  concluded  that  departures  of  bright¬ 
ness  or  radiance  in  small  areas  from  the  average  value  within  the 
cloud  cover  over  a  large  area  could  be  used  to  infer  precipitation. 
Cheng  and  Rodenhuis  (1977)  used  NOAA  2  (National  Oceanic  and 
Atmospheric  Administration  2)  visible  and  infrared  imagery  along  with 
WSR-57  radar  data  in  comparing  satellite  images  with  radar  rainfall 
rates.  They  found  correlation  coefficients  between  satellite  data 
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and  radar  rainfall  rates  of  from  0.23  to  0.34  for  visible  data  and 
from  0.17  to  0.27  for  infrared  data. 

Sil'dar  (1972)  evaluated  the  relationship  of  rainfall  patterns  to 
cloud  brightness  fields  and  found  that  all  of  the  radar  echoes  w-.  re 
included  in  the  upper  3G  to  40  per  cent  of  the  ATS  3  brightness 
range  for  the  two  cases  considered.  He  indicated  an  excellent  cor¬ 
relation  between  positive  brightness  anomaly  area,  radar  echo  area, 
and  the  precipitation  pattern.  Reynolds  and  Vender  Ilaar  (1973) 
claimed  that  when  growth  of  clouds  could  be  detected  by  radar  a 
corresponding  increase  in  reflected  radiance  could  be  detected  in 
ATS  3  visible  imagery.  Griffith  and  Woodley  (1  973)  compared  re¬ 
flected  radiances  with  radar-measured  echo  heights  and  found  a 
positive  correlation  between  cloud  brightness  and  echo  height. 
However,  Negri  et  al.  (1970)  reported  that  brightness  did  not  seem, 
to  correlate  well  with  any  echo  parameters.  Negri  et  a  1 .  also  indi¬ 
cated  that  cloud  area  changes  led  echo  area  changes  for  the  storm 
they  examined.  Adler  and  l?onn  (1977)  found  that  duration  of  cloud 
elements,  minimum  cloud-top  temperature,  and  cold  area  rate  of 
increase  seemed  to  be  correlated  with  the  occurrence  of  severe 
weather. 

Reynolds  (1978)  used  an  interactive  computer  system  to  evaluate 
digital  radar  and  digital  GOES  data.  Reynolds  indicated  that  infrared 


ckita  show  relative  changes  in  "iop  height  which  are  related  well  to 
rainfall  and  that  maximum  rainfall  occurs  before  the  coldest  cloud- 
top  temperature  is  reached.  Reynolds  ot  al.  (1978)  also  considered 
the  use  of  meteorological  satellites  in  weather  modification  projects 
and  included  comparisons  of  radar  data  and  geostationary  satellite 
data  in  seeking  further  to  explore  satellite-rainfall  relationships. 
Smith  and  Reynolds  (1978)  have  developed  procedures  to  produce 
numerical  composites  of  digital  radar  and  digital  satellite  data  for 
use  in  analyzing  precipitation  patterns. 

Griffith  et  aj_.  (1978)  employed  ATS  3  and  SMS/GOES  imagery  in 
rainfall  estimation  applications  using  digital  ladar  and  roinqage  data 
for  ground  truth.  They  developed  a  rainfall  estimation  pro¬ 
cedure  based  on  the  growth  of  radar  echoes  and  satellite  cloud  areas 
and  examined  the  procedure  for  potential  real-time  applications  and 
automation . 

3.  Objectives  of  the  Investigation 

The  primary  objective  of  this  investigation  is  to  compare  digital 
meteorological  radar  data  with  infrared  imagery  from  a  meteorological 
satellite  in  geostationary  earth  orbit  by  determining  the  correlation 
between  point  values  of  infrared  digital  count  and  point  values  of 
radar  reflectivity  and  vertically  integrated  liquid  water  content  (VIL). 
Comparisons  are  made  between  infrared  digital  count  from  the 


satellite  sensor  and  zero-tilt  reflectivity,  three  levels  of  partial 
vertically-summed  reflectivity  (PVSZ) ,  and  VI L,  for  the  purpose  cl 
determining  which  sets  of  data  demonstrate  the  best  relationship 
between  the  satellite  and  radar  data. 

A  secondary  objective  of  this  investigation  is  to  develop  and 
implement  procedures  for  computer  processing  of  GOES  imagery  stored 
on  magnetic  tape  and  to  map  the  data  to  points  in  the  Texas  A & M 
University  (TAMU)  digital  weather  radar  grid. 

4.  Techniques  and  Scope  of  the  Investigation 

The  radar  data  for  this  investigation  were  obtained  from  the 
digital  weather  radar  system  in  the  Department  of  Meteorology  at 
TAMU.  The  reduction,  interpolation  and  display  of  the  radar  data 
were  accomplished  on  the  Amdahl  470  V/6  computer  at  TAMU  using 
the  technique  developed  by  Sieland  (1977)  as  modified  by  McAnelly 
(1979).  The  digital  GOES  infrared  data  were  obtained  from  the  Man- 
computer  Interactive  Data  Access  System  (McIDAS)  at  the  Space 
Science  and  Engineering  Center  of  the  University  of  Wisconsin  in 
Madison,  Wisconsin.  Reduction  and  display  of  the  digital  satellite 
data  also  were  accomplished  by  the  Amdahl  470  V/G  computer  using  a 
technique  designed  to  map  the  satellite  data  to  individual  points  in 
the  radar  grid  established  by  the  Sieland  technique. 


Tilt-sequence  10-orri  digital  radar  data  were  collected  iv-ar- 
simultaneously  with  GOES  infrared  imagery  on  3  May  197S.  The  tilt- 
sequence  data  were  processed  to  produce  point  values  of  zero-tilt 
reflectivity,  three  levels  of  partial  vertically  summed  reflectivity 
(PVSZ) ,  and  vertically  integrated  liquid  water  content  f'/JL).  The 
point  values  of  the  digital  radar  data  were-  then  tabulated  along  with 
GOES  infrared  digital  count  data  mapped  to  corresponding  points. 

The  GOES  infrared  data  were  mapped  from  images  at  times  correspond¬ 
ing  to  the  start  of  the  radar  tilt-sequences  and  the  end  of  the  tilt- 
sequences.  The  tabulated  data  were  punched  onto  computer  cards 
and  processed  to  compute  correlations  between  each  field  of  radar 
and  satellite  data. 

Three  sets  of  radar  and  satellite  data  were  processed  to  reduce 
computer  costs  and  meet  time  constraints  imposed  by  tabulating  the 
data  and  punching  the  data  onto  computer  cards.  One  set  of  data 
was  selected  during  the  early  development  of  radar  echoes  in  the 
TAMU  area  and  included  one  tilt-sequence  of  radar  data  and  two 
infrared  images  for  a  corresponding  area.  The  remaining  two  data 
sets  were  selected  from  a  time  when  radar  echo  activity  in  the  area 
appeared  to  be  at  or  near  its  maximum  intensity.  The  last  two  data 
sets  were  taken  from  a  single  tilt-sequence  and  from  two  correspond¬ 
ing  infrared  images. 


CHAPTER  II 


METEOROLOGICAL  RADAR  DATA 

Microwave  radar  is  one  ol  the  many  tools  employed  by  meteoro¬ 
logists  in  studying  the  atmosphere.  Digital  radar  systems  have  made 
possible  quantitative  companions  of  large  volume:  f  radar  data  with 
data  from  other  sources  including  the  mi  t. •> olegk;  !  satellite  data 
employed  in  this  study.  Much  of  the  mV  'L  l  in  this  chapter 
concerned  with  basic  radar  theory  and  th-  TAMU  weather  radar  system 
is  extracted  from  Neyland  (1  978)  that  contains  <5  me  re  complete 
description  of  the  digital  radar  system  and  f  ho  Sic  Mad  (1377)  com~ 
puter  program. 

1.  Basic  Radar  Theory 

The  basic  form  of  the  radar  equation  used  by  radar  meteorologists 
was  derived  by  Probert- Jones  (1962).  Probcrt-Jones  assumed  a  more 
realistic  beam  shape  than  had  been  used  in  previous  derivations  and 
further  assumed  that  no  significant  attenuation  existed  between  the 
radar  antenna  and  the  target  and  that  the  target  completely  filled  the 
spatial  volume  illuminated  by  the  radar  beam.  With  these  assump¬ 
tions,  the  average  backscattered  power,  Pr  (watts),  received  from 
the  target  at  range  r  (km) ,  is  given  by 
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F>  =  cJ’rJ  —  z. 


(l) 


where  c  (w  km^  in^  mm  is  the  radar  constant  unique  to  the  indivi¬ 
dual  radar,  |kp  (dimensionless)  is  the  dielectric  constant  from 
scattering  theory,  and  ZG  (mmD  m  °)  is  the  equivalent  radar  reflec¬ 
tivity  factor.  An  expression  for  the  radar  constant  is  given  by 


c  =  fr3  Pt  G2  h6  (j)  (2) 

512  A2  2  In  2 

where  Pj.  ;.s  the  transmitted  power  during  a  radar  pulse  (w),  G  is  the 
antenna  gain  (dimensionless),  h  is  the  pulse  length  in  space  of  the 
transmitted  pulse  (cm),  0  is  the  horizontal  beamwidih  (radians),  (J) 
is  the  vertical  beamwidth  (radians),  and  A  is  the  wavelength  of  the. 
transmitted  microwave  energy  (cm).  The  factor  2  In  2  is  the  Frobert- 
Joncs  correction  factor.  The  value  of  c  is  8.G09  x  10_^  for  the 
TAMU  10  cm  radar  and  the  value  of  |k|3j  is  0.93  (Battan,  1  973). 

The  output  of  the  radar  system  is  digitized  values  of  returned 
power,  Pr  ,  which  must  be  converted  to  values  of  reflectivity  before 
further  data  reduction.  The  usual  practice  is  to  measure  the  returned 
power  in  terms  of  decibels  with  respect  to  a  standard  reference  power 
level,  normally  1  mw.  Power  levels  are  then  expressed  in  units  of 
dBm,  cither  above  (+)  or  below  (-)  1  mw  according  to 


P(dBm)  =  10  log  ,  n  . 

1  g-3(watts) 


(3) 
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To  determine  values  of  Ze,  Eq.  (1)  is  solved  for  Ze  giving 


Z  =  -I— 

6  c  |k)2 


Taking  the  logarithm  of  Eq.  (4)  gives 


(4) 


logiQ  (Ze)  =  2  log1Qr  +  log]Q  Pr  -  log10c  J k| 2  .  (5) 

For  the  TAMU  1  0  cm  radar,  c  is  equal  to  8 . 609  x  1  0~^  ^  and 
logiocIM“  equals -10.1.  Substituting  these  values  into  Eq.  (5)  gives 

log10Ze  =  2  log10r  +  logio  pr  +  10  - 1  •  (6) 

The  digital  value  of  P  is  converted  to  its  dBm  equivalent  (always 
negative)  through  the  use  of  calibration  data.  The  received  power 
in  watts  is  related  to  the  received  power  in  dBm  by 

log10Pr  (watts)  =  0.1  Pr  (dBm)  -  3,0.  (7) 

Substitution  of  Eq.  (7)  into  Eq.  (6)  gives 

logi  q  ZQ  =  0.1  Pr  +  2  log1 0  r  +  7 . 1 .  (8) 

Finally,  the  value  of  Zc  is  given  by 

Z0  -  10(°'l  Pr  +  21og10r  +  wc)  (9) 

where  Ze  is  in  units  of  mm^  m-2  and  wc  is  7.1 . 

The  several  assumptions  inherent  in  the  development  of  Eq.  (9) 
include: 

1)  The  transmitted  microwave  energy  is  not  significantly 
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attenuated  between  the  radar  antenna  and  the  target.  This  assump¬ 
tion  has  been  shown  to  be  valid  for  radars  with  wavelengths  near 
10  cm  (Greene,  19G4). 

2)  The  spatial  volume  illuminated  by  the  radar  beam  is  com¬ 
pletely  filled  by  the  target.  This  condition  is  not  always  met, 
particularly  on  the  periphery  of  a  storm,  with  a  resulting  loss  of 
resolution  of  fine  scale  detail  in  such  areas  (Greene,  1971).  How¬ 
ever,  this  effect  does  not  significantly  affect  the  major  features  in 
the  digital  data . 

3)  The  Rayleigh  approximation  is  used  to  describe  the 
scattering  properties  of  spherical  liquid  water  drops  having  diameters 
on  the  order  of  0.04  A,  where  A  is  the  radar  wavelength  in  cm.  In 
severe  storms,  the  large,  non-sphcrical  water  drops  and  hailstones 
which  may  be  present  do  not  meet  the  Rayleigh  criteria.  However, 
the  enhanced  reflectivity  of  those  particles  may  be  useful  in  identi¬ 
fying  such  storms. 

4)  Each  digital  datum  value  represents  a  point  in  the  center 
of  the  radar  volume. 

5)  The  equivalent  radar  reflectivity  values,  Z0,  obtained 
from  the  digital  data,  are  representotive  of  a  continuous  scalar  field. 

2 .  Earth  Curvature  Correction 

Microwave  radiation  propagating  in  free  space  will  follow  a 


straightline  path.  However,  the  microwave  radiation  of  a  radar  beam 
in  the  atmosphere  does  not  follow  a  straightline  path  but  instead 
follows  a  curved  path  due  to  refraction.  The  amount  of  curvature  of 
the  path  depends  on  the  index  of  refraction  (n).  Under  normal 
atmospheric  conditions  where  temperature  and  humidity  decrease 
with  height  the  path  of  the  radar  beam  curves  slightly  downward  to¬ 
ward  the  surface  of  the  earth. 

Ray  theory  may  be  applied  to  the  problem  if  the  change  with 

height  of  the  index  of  refraction,  — — ,  is  small.  For  that  case,  the 

dh 

exact  differential  equation  for  a  ray  in  a  spherically  stratified  at- 
rnosph  is 


d2h  _  (  2  1  dn  )  (dli) 2  +  ,R  + 1) ,  2  .  1  +  1  dn  >  (10) 

ds2  R  +  h  n  dh  ds  R  R  +  h  n  dh 


whore  h  is  the  height  of  the  beam  above  the  earth’s  surface  at  a 
distance  s  from  the  transmitter,  R  is  the  earth's  radius,  and  n  is  the 
index  of  refraction.  Since  £,  the  angle  at  which  the  beam  is  sent 


out  measured  from  a  horizontal  plane,  is  usually  very  small, 

dh  2  o 

(— r— )  =  tan  £  Z.  /II .  Also,  since  n  ~  1  and  hz^ZR,  Eq.  (1 0)  can  be 
ds  ~ 

reduced  to 


d2h  _  J_  + 
ds2  R  dh 


(ID 


In  this  investigation  it  is  necessary  to  consider  the  radar  beam  axis 


as  a  straight  lino  so  that  a  fictitious  earth  radius  given  by 


1  +  dn 


must  be  assumed  (Rattan,  1  973).  Since  is  smell  and  nearly  linear 

dh 

with  a  value  of  -4  x  10~^  m-^  ,  an  earth  curvature  correction  of 


R‘  =  ~  R  is  used  and  the  resultant  radar  beam  axis  may  be  considered 
3 

to  be  a  straight  line. 


3.  The  TAMU  Weather  Radar  System 

The  TAMU  Weather  Radar  System  consists  of  both  a  3-cm  and  a 
10-cin  radar  operating  in  parallel.  Each  of  the  two  radar  systems 
has  a  nearly  identical  makeup  and  consists  of  the  following  major 
sub-systems:  antenna,  conventional  analog  weather  radar  set,  di¬ 
gital  video  integrator  and  processor  (DVIR)  with  plan  position  indica¬ 
tor  (PPI)  display,  and  a  nine-track  tape  unit  to  record  the  buffered 
digital  output  of  the  DVIP.  Since  only  10-cm  radar  data  are  used  in 
this  study,  the  following  discussion  is  limited  to  some  of  the 
characteristics  of  the  10-cm  radar. 

The  10-cm  radar  was  constructed  at  TAMU  using  components  from 
several  different  radars.  The  radar  receiver  is  mounted  in  a  CPS-9 
console  and  uses  the  same  synchronizing  pulse  generator  used  to 
trigger  the  3-cm  radar.  This,  combined  with  the  common  mounting 
of  the  3- and  lo-cm  antennas,  permits  simultaneous  scanning  of  the 
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same  illuminated  volume  by  both  radars. 

The  DVIP  is  a  high  speed  data  acquisition,  digitizing,  and  pro¬ 
cessing  system  which  continuously  averages  radar  logarithmic  video 
output  in  range  (range  averaging)  and  in  direction  of  antenna  scon 
(time  averaging)  using  exponentially  weighted  digital  integrator 
techniques  and  synchronized  by  the  radar  system  trigger.  The  DVIP 
installation  provides  for  real-time  contoured  PPI  displays  and  digital 
recording  of  the  digital,  integrated  data. 

The  DVIP  output  signal  is  split  into  two  channels,  the  digital 
channel  and  the  display  channel.  The  display  channel  provides 
either  log  video  or  contoured  log  video  to  the  PPI  scope.  The  digital 
output  consists  of  the  1-  or  2-km  range  increment  integrated  video 
samples,  provided  in  a  buffered  8-bit  parallel  binary  configuration. 
Each  digital  output  word  represents  a  digital  value  of  the  integrated 
video  intensity  on  a  linear  scale  over  the  input  dynamic  range  of 
80  dB.  At  the  digital  output  there  are  215  2-km  or  430  1-km  digital 
samples  representing  the  range  increments  from  21  through  450  km 
of  range,  depending  on  the  range  increment  selected. 

The  digital  output  channel  of  the  DVIP  is  fed  into  a  9-track 

0 

magnetic  tape  recording  unit.  This  tape  drive  unit  records  t he  digital 
data  in  a  format  that  makes  it  readily  accessible  by  computer  and 


controls  the  azimuth  increment  integration  of  the  data. 


4.  Data  Reduction  and  Display 

The  digital  output  of  the  DVIP  is  recorded  on  9 -track  magnetic 
tape;  which  becomes  the  input  to  the  TAMU  Amdahl  470  V/6  computer. 
Computer  processing  of  the  digital  radar  data  is  performed  using  a 
program  developed  by  Sieland  (1977)  and  modified  by  IVcAnelly 
(1979). 

The  computer  program  begins  by  transforming  the  spherical  co¬ 
ordinate  system  (r,Cx,£)  of  the  recorded  digital  radar  data  into  a 
cylindrical  coordinate  system  h)  with  three  separate  height 

divisions,  or  partial  vertically  summed  reflectivity  (PVSZ)  layers. 

The  data  are  then  transformed  into  a  two-dimensional  rectangular 
coordinate  system  (x,  y)  and  five  maps  are  produced,  a  zero  level 
(zero-tilt)  reflectivity  map,  the  three  PVSZ  maps,  and  a  map  of 
vertically-integrated  liquid  water  content  (VIL) .  The  conversion  of 
coordinate  systems  requires  tire  use  of  two  interpolation  schemes, 
a  Lagrungian  linear  or  cubic  interpolation  scheme  along  each  radial 
of  data  and  a  quadratic  scheme  for  the  conversion  to  rectangular 
coordinates.  An  extensive  explanation  of  these  interpolation  schemes 
is  contained  in  Sieland  (1977). 

The  three  PVSZ  maps  arc  produced  by  summing  the  reflectivity 
values  in  a  vertical  column  above  a  given  (x,  y)  point  at  the  surface 
into  the  lower  layer  (0-deg  to  15  kft),  the  middle  layer  (15  to  35 


18 


kft),  and  the  upper  layer  (38  to  50  kit).  The  .summed  values  of  re¬ 
flectivity  in  each  layer  at  each  (x,  y)  point  are  divided  by  the  number 
of  tilt  angles  used  to  produce  the  sum  so  that  a  mean  value  of  the 
reflectivity  is  produced,  VIL  is  produced  using analgorithm  developed 
by  Greene  (1971)  and  is  based  on  the  relationship 

i~  r  htop  4/7 

M*  =  3.44  >:  10~6  \  ze  dh  (13) 

-J  h, 

bottom 

O 

wlieie  M*  is  the  VIL  in  units  of  kg  m  and  b  is  the  height  in  meters. 

Because  of  the’  large  amount  of  data  processed  and  to  reduce 
computation  time,  the  Sicland  program  maps  only  a  selected  portion 
of  a  grid  centered  on  TAMU.  The  mapped  portion  is  selected  by 
specifying  a  distance  east  (-)  or  west  (+)  from  TAMU  to  the  lower 
left,  corner  of  the  100  x  100  km  grid,  a  distance  north  (--)  or  south  (+) 
from  TAMU  to  the  lower  kft  conv  of  the  100  x  100  kin  grid,  a 
starting  azimuth  for  data  processing,  and  a  last  azimuth.  Computer 
processing  of  the  data  on  the  100  x  100  km  grid  results  in  a  data 
output  spatial  resolution  of  2-km  on  a  51  x  51  output  array. 

The  radar  data  used  in  this  investigation  were  produced  and  re¬ 
corded  by  the  TAMU  Weather  Radar  System  during  the  evening  of  2 
May  197  8  (3  May  1978  GMT),  as  an  area  of  light  and  moderate 
thunderstorm  activity  moved  through  the  College  Station  area.  Tilt- 
sequence  data  were  recorded  at  times  which  would  permit  comparison 


with  GOES  satellite  data  colleeU  d  that  evening  at  IS -min  interva ; v . 
This  investigation  used  two  till. --.sequences  of  the  radar  data  for 
comparison  with  COES  imagery,  the  first  tilr-sequenco  beginning  at 
0022  GMT  and  the  tenth  tilt-sequence  beginning  at  0 2 1  E  GMT.  Tin  sc 
tilt-sequences  were  processed  by  Mr.  Ray  McAnelly  to  produce  maps 
of  zero-tilt  reflectivity,  PVSZ  in  three  layers,  and  VI L. 

The  maps  of  zero-tilt  reflectivity,  PVSZ,  and  7IL  wore  adjusted 
for  ground  clutter  effects  by  excluding  all  data  within  a  -10  km  radius 
of  TAMU  in  order  to  avoid  data  contamination  and  to  assure  compara¬ 
bility  between  results  for  the  r.eio-tilt  reflectivity  maps  and  the  PVSZ 
and  VII,  maps.  Data  at  the  outer  edges  of  the  radar  maps  end  along 
bounding  azimuths  used  by  the  Sic  land  program  also  were  excluded 
duo  to  suspected  boundary  into  polation  problems  in  the  computer 
program.  After  these  adjustments,  the  da  Us  were  tabulated  for 
comparison  with  the  infrared  satellite  data. 
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CHATTER  III 

GEOSTATIONARY  METEOROLOGICAL  SATELLITE  DATA 

Meteorological  satellites  have  boon  used  by  meteorologist;.:  for 
both  qualitative  and  quantitative  investigation?:  of  atmospheric 
phenomena.  Geostationary  satellites,  by  virtue  of  near-continuous 
coverage  of  large  areas  of  the  earth  surface,  permit  more  detailed 
examination  of  severe  weather  dev<  lopmont  than  :'s  possible  through 
the  use  of  polar-orbiting  satellites.  The  on-board  digitisation  of  the 
data  from  the  newer  geostationary  satellites  makes  possible.-  ready 
comparison  of  those  data  with  other  digital  data  fields  and  is  usefuj  ir 
the  analysis  of  severe  weather  occurrences. 

1.  The  S MR/GOES  Spacecraft 

The  Synchronous  Meteorological  Satellite/ Geostationary 
Operational  Environmental  Satellite  (SMS/GOES)  spacecraft  is  a 
cylindrical  spacecraft  weighing  approximately  5  78  kg  at  launch. 

The  spacecraft  arc  placed  into  a  circular  orbit  at  about  35,800  km 
above  the  earth  surface.  The  satellites  orbit  the-  earth  at  a  speed 
c  about  1  1,000  km  hr  .  This  combination  of  altitude  and  speed 
permits  each  satellite  to  remain  continuously  above  the  same  point 
on  the  earth  and  produces  what  is  referred  to  as  a  geostationary  or 
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geosynchronous  obit . 

The  spacecraft  is  controlled  for  proper  earth-imaging  by  an 
attitude  control  subsystem  which  maintains  the  spin  rate  at  ICO  rpm 
and  aligns  the  spacecraft  spin  axis  parallel  to  the  earth's  polar 
axis.  On-board  altitude  sensors  are  used  to  determine  the  spin 
axis  orientation  and  a  hydrazine  jet  system  is  used  to  adjust  the 
spin  axis  orientation,  spin  rate,  and  orbital  position  as  necessary 
(Corbell  et  nl_.  ,  197C). 

The  primary  purpose  of  trie  spacecraft  is  to  provide  earth-imaging 
in  the  visible  spectrum  (0.55  to  0.7p.m)  with  an  0.9-kin  resolution  at 
subpoint  and  in  the  infrared  spectrum  (10.5  to  12.6p.nn)  with  an  8-kn 
resolution  at  subpoint.  Tire  spacecraft  also  possesses  a  capability 
to  provide  direct,  quantitative  measurements  of  solar  activity  and 
to  collect  and  distribute  environ  mental  data  measured  on  remotely 
located  data  collection  platforms  on  the  earth's  surface  and  in  its 
atmosphere . 

Primary  sub-systems  aboard  the  spacecraft  include  the  Visible 
and  Infrared  Spin  Scan  Radiometer  (VISSR),  the  Space  Environmental 
Monitor  (SEM),  the  Telemetry,  Tracking,  and  Command  (TTC)  sub¬ 
system,  and  an  attitude-control  sub-system.  The  VISSR  provides 
earth-imaging  in  the  visible  and  infrared  spcctrums.  The  SEM  in¬ 


cludes  a  magnetometer,  a  solar  x-roy  telescope,  and  an  energetic 
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particle  monitor.  The  TTC  includes  equipment  for  S-bnnd  tj. ins- 
mission  and  reception  (in  reduced  bandwidth)  of  VISSR  data,  S-band 
transmission  of  weather  facsimile  data,  UI1F  reception  of  da  la  from 
remote  data  collection  platforms  and  downlink  of  these  data  to  earth, 
and  VI IF  and  S-band  communications  equipment  for  commanding  the 
spacecraft,  for  telemetry,  and  for  transmitting  SUM  data.  The 
primary  instrument  providing  the  data  used  by  most  meteorologists 
is  the  VISSR  (Corbcll  ot  ab  ,  197G). 

2.  The  Visible  and  Infrared  Spin  Scan  Radiometer 

The  Visible  and  Infrared  Spin  Scan  Radiom.  ter  scans  from  west 
to  east  in  eight  identical  visible  channels  and  two  infrared  channels. 
The  sensor  provides  visible  data  at  0.9 -km  resolution  at  satellite 
subpoint  and  infrared  data  of  8- km  resolution  ot  satellite  subpoint. 
The  resolution  in  both  the  visible  and  intruded  data  deteriorates  for 
points  away  from  the  satellite  subpoint.  With  the  satellite  rotating 
at  100  rpm,  the  VISSR  scans  the  earth  for  about  one-twentieth  of  each 
complete  rotation.  The  radiometer  performs  1821  steps  in  successive 
scans  from  north  to  south  in  10.2  min  and  provides  an  image  of 
the  complete  earth  disc  or  about  one  quarter  of  the  earth's  surface. 
The  resulting  visible  images  contain  14,5  08  lines.  The  infrared 


images  have  a  total  of  1821  lines.  In  addition  to  the  normal  scanning 
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mode.,  the  spaced  aft  may  be  placed  into  a  limited  scan  mode .  hi 
the  limited  scan  mode  the  north-south  soon  may  be  limited  to  a 
fewer  number  of  scan  lines,  thereby  reducing  the  area  of  coverage 
but  increasing  the  frequency  of  imaging. 

The  VISSR  senses  radiation  in  both  the  visible  and  infrared 
spoclrums.  Radiation  is  received  by  the  VISSR' s  primary  optics  via 
a  45°  object-space  scan  mirror.  The  mirror  is  an  eiliptically  shaped 
plane  mirror  which  is  tilted  about  its  minor  axis  to  obtain  the  north- 
south  scan  stops.  West-east  scans  are  produced  by  spacecraft 
rotation  about  the  spin  axis.  Energy  from  the  scan  minor  is  collected 
by  a  Ritchey-Clu alien  optical  system  which  includes  a  baffle  that 
extends  from  the:  mirror  to  minimize  the  effects  of  scattered  radiation . 
Radiation  in  the  visible  spectrum  is  detected  at  the  prime  focus  using 
eight  fiber  optics  in  a  linear  array  at  the  focal  plane.  The  oilier  ends 
of  the  fibers  are  optically  integrated  with  eight  photomultiplier  tubes 
having  the  desired  0 . 55- to  0.70-jum  response.  The  prime  focal 
plane  is  also  relayed  to  the  long-wavelength  HgCd/Te  detector  using 
two  germanium  relay  lenses.  An  optical  filter  between  the  final  re¬ 
lay  lens  and  the  detector  restricts  the  radiation  to  the  1 0 . 5-  to  1 2 . 6- 
pm  wavelength  band  (McKowan,  1977). 

3.  Data  Distribution  and  Processing 

The  SMS/GOES  VISSR  data  are  distributed  as  indicated  in  Fig. 


1.  Row  data  arc  rood  out  and  prcprocessod  at  the  Command  and  Data 
Acquisition  (CDA)  station  at  Wallops  Island,  Virginia.  The  CDA 
processes  the  data  to  reduce  the  data  rate  hv  about  16  to  1  for  simpli¬ 
fication  of  data  transmission  and  then  transmits  those  stretched  data 
bock  to  the  satellite.  Lower  resolution  infrared  data  are  formatted 
for  analog  transmission  and  are  sent  via  telephone  lines  from  the 
CDA  to  the  Satellite  Field  Service  Stations  (SPSS)  and  to  the  Central 
Data  Distribution  Facility  (CDDF)  at  Marlow  Heights,  Maryland. 

The  stretched  data  are  retransmitted  by  the  satellite  to  the 
National  Environmental  Satellite  Service  (NESS)  facility  at  Suitlund, 
Maryland, and  to  other  facilities  with  a  direct  readout  capability. 

The  NESS-Suitiand  complex  relays  the  stretched  data  via  microwave 
link  to  the  CDDF.  The  CDDF  formats  the  data  and  prepares  it  for 
transmission  to  SPSS's  throughout  the  country.  Enhancement  of  the 
infrared  data  to  emphasize  specific  features  in  the  imagery  can  be 
accomplished  both  by  the  CDA  at  Wallops  Island  and  by  the  CDDF. 

GOES/SMS  VISSR  data  are  available  to  many  users  through  a 
"GOES-Tap"  program  using  high-quality  telephone  line  transmission 
of  the  data  from  an  SFSS  to  the  user.  As  an  alternative,  some  faci¬ 
lities  have  developed  a  direct  readout  capability  for  receiving  the 
stretched  VISSR  data  (Corbell  ct  aj_.  ,  1976). 


4.  Basic  Characteristics  of  Infrared  Satellite  Data 

Infrared  radiation  is  omitted  by  objects  at  an  intensity  that  de¬ 
pends  on  the  temperature  of  the  body.  If  the  object  is  assumed  to 
emit  radiation  as  a  "black  body,"  me.,  to  absorb  all  radiation  falling 
upon  it  without  reflecting  any  of  the  incident  radiation,  then  the 
monochromatic  intensity  of  the  emitted  black  body  radiation  is  given 
by  Planck's  Law, 

E  =  3.74  x  1 0“5  A~5  (exp  — -l)'1  (1) 

where  A  (cm)  is  the  wavelength  of  the  emitted  radiation  and  T  (K  ) 
is  the  temperature  of  the  black  body  (Ilaltiner  and  Martin,  1957). 

The  proportionality  between  the  temperature  of  the  target  and  the 
intensity  of  the  infrared  radiation  emitted  by  the  target  forms  the 
basis  for  all  infrared  imaging. 

Infrared  sensors  aboard  meteorological  satellites  measure  long¬ 
wave  (usually  1  0 . 5- to  1 2 .G-yjm  wavelength)  radiation  emitted  by 
cloud,  land,  and  water  surfaces.  These  measurements  may  be  con¬ 
verted  to  black  body  temperatures  which  generally  represent  the 
temperatures  of  the  surfaces  viewed.  This  characteristic  of  infrared 
satellite  imagery  is  especially  useful  in  that  it  permits  an  approxi¬ 
mate  evaluation  of  the  height  of  various  clouds  when  the  atmospheric 
temperature  is  known  from  upper-air  data.  The  association  of  cold 


cloud  tops  apparent  in  infrared  da!-),  convective  activity,  and  pre¬ 
cipitation  suggested  by  some  of  the  authors  cited  in  Chapter  I  forms 
the  basis  for  much  of  the  following  material  in  this  investigation. 

5.  Mel  DoS  Data 

The  Mon -computer  Interactive  Data  Access  System  (Mel DP  S) 
was  developed  by  the  Space  Science  and  Engineering  Center  (SSEC) 
at  the  University  of  Wisconsin  os  a  means  of  obtaining  easy  access 
to  satellite  data  in  the  time  domain  for  examining  atmospheric  wea¬ 
ther  systems  as  dynamic  systems  (Smith,  1975).  The  operator  using 
Mol  DAS  views  image  and  graphics  data  on  a  large  color  television 
monitor.  Images  may  be  viewed  in  single  trainee  or  sequences  and. 
may  be  enhanced,  colori ;v  <1 ,  or  combined  at  the  choice  of  the 
operator.  Control  of  the  system  is  effected  from  the  operator's  con¬ 
sole  by  means  of  twin  joysticks  nd  a  keyboaid.  The  operator  also 
may  be  supported  by  hard  copy  generating  devices  such  as  line 
printers  and  graphics  plotters  (Haig,  1  978). 

Image  data  used  by  the  McIDAS  also  may  be  saved  on  magnetic 
tape  for  later  use  or  foi  the  use  of  oilier  investigators.  The  imagery 
used  in  this  study  was  produced  in  such  a  fashion  by  the  McIDAS  at 
the  SSEC  in  Madison,  Wisconsin. 

McIDAS  image  save  tapes  contain  two-dimensional  image  sectors 
prepared  from  a  variety  of  image  sources.  The  data  for  this  study 
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arc  a  sequence  of  infrared  images  from  the  GOFS-East  satellite 
taken  between  0015  GMT  and  0345  GMT  on  3  May  1978.  The  data 
are  stored  on  the  tapes  in  a  directory  record  which  contains  date, 
time,  and  line  and  element  coordinate  information  necessary  for 
image  processing  and  in  image  sector  records  which  follow  the  di¬ 
rectory  record.  Each  image  sector  record  contains  six  lines  of 
GOES  data  in  24-bit  words  so  that  for  a  complete  image  of  5 Of.'  lines 
a  total  of  84  image  sector  records  is  required.  The  directory  record 
is  225  words  long  and  each  image  save  record  is  2211  words  long 
(G723  8-bit  bytes).  6720  bytes  of  each  image  save  record  are  used 
for  data  with  the  remaining  3  bytes  containing  a  lceoul  sequence- 
counter  (Young,  1978). 

Full-resolution  GOES-East  infrared  imagery  are  used  in  this 
study.  Each  pixel,  or  picture  element,  has  a  resolution  of  approxi- 
•mately  8  km  in  the  line  direction  (N-S  at  satellite  subpoint)  and  4  km 
in  the  element  direction  (E-W  at  satellite  subpoint).  The  Mol  DAS 
artifically  subdivides  each  pixel  to  a  4  km  by  4  kin  pixel  on  the 
image  save  tapes.  Line  and  element  numbering  in  the  infrared  data 
is  based  on  the  system  used  for  full-resolution  visible  data 


(Young,  1978). 


6.  Data  Reduction  and  Display 


Use  of  the  McIDAS  image  save  tapes  requited  development  of 
procedures  to  navigate,  or  locate,  the  image  pixels,  to  map  the 
infrared  data  into  the  grid  used  for  the  digital  >vidar  data,  and  to 
correct  the  data  for  displacement  error  due  to  cloud  height, 

a.  Imago  Navigation 

Navigating  the  image  pixels  requires  extensive  computer  soft¬ 
ware  that  was  not  available  for  this  investigation.  Instead,  Mr. 

J.  T.  Young  of  the  SSEC  provided  line  and  element  navigation  solu¬ 
tions  for  four  latitude  and  longitude  points  in  a  grid  around TAMU 
for  each  hour  from  0000  GMT  to  0400  GMT  on  3  May  1978.  A  fourth- 
order  polynomial  was  than  fitted  to  the  line  and  element  solutions 
for  eacii  point  to  permit  interpolating  line  and  element  solutions  for 
intermediate  times.  The  interpolated  line  and  element  solutions 
then  were  used  as  inputs  to  a  computer  program  called  MAC  MAP  3 
which  produced  mapping  equations  to  map  the  satellite  data  into  the 
radar  grid . 

b.  Mapping  the  Satellite  Data  to  the  Radar  Grid 

The  satellite  data  array  is  a  two-dimensional  array  arranged  by 
line  and  element  number.  The  radar  grid  ccnteicd  on  TAMU  can  be 
regarded  as  an  (x,  y)  array  with  a  positive  y-oxis  oriented  due 
North  and  a  positive  x-axis  oriented  due  Easl.  The  variation  in 


line  find  element  solutions  lor  the  four  points  navigated  by  the  SSCC 
suggested  that  simple  bi- linear  relationships  of  the  form 

line  number  -■  a  x  -i  b  y  +  c  (2) 

element  number  -  dx  +  ey  +  f  (3) 

could  be  used  to  relate  the  digital  satellite  line  and  element  numbers 
to  the  coordinates  of  the  radar  gi id, where  the  constants  a,  b,  c,  d, 
e,  and  f  arc  determined  by  solving  three  simultaneous  equations  using 
solutions  for  three  cf  the  four  points  navigated  by  the  SSCC.  The 
original  mapping  program,  MAC  MAP  1,  used  navigation  solutions  for 
throe  points  to  produce  coefficients  for  the  mapping  equations  and 
used  the  remaining  solution  for  the  fourth  point  to  choc),  the  accuracy 
of  the  equation.  MACMAP  1  solutions  produced  results  accurate  to 
one  half  pixel  for  the  fourth  point  for  all  images  except  two  and  was 
within  one  pixel  for  those  two  images. 

An  effort  to  produce  more  accurate  mapping  of  the  satellite  data 
resulted  in  adding  new  terms  to  the  mapping  equations  using  the 
fourth  point  to  determine  "correction  factors"  and  );2.  The  new 
equations  wore  of  the  form 

lino  number  =  ax+by+cH  k^(x-x0)(y-y0)  (4) 

element  number  =  dx  +  ey  +  f  +  k2  (x-xQ)  (y-y0)  (5) 

where  the  coefficients  a,  b,  c,  d,  c,  and  f  arc  determined  by  solving 
simultaneous  equations  for  the  three  original  points  and  the  correction 
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factors  by  determining  the  difference  between  the  SSI.'C  navigation 
solutions  and  the  solutions  from  Fq.  (2)  ami  Fq.  (3).  The  values  of 
the  constants  >:Q  and  yQ  ore  determined  from  the  "common  values" 
of  x  and  y  shared  by  the  three  points  used  to  solve  the  simple  bi¬ 
linear  equations.  Fig.  2  illustrates  the  gcneia!  arrangement  of  the 
four  navigated  points  relative  to  TAhiU  and  the  radar  grid  and  illu¬ 
strates  how  x0  and  yQ  were  determined. 

The  MACMAP  3  computer  program  produced  coefficients  foi  fq. 

(4)  and  Fq.  (5)  for  each  of  the  images  used  in  this  study.  The 
solutions  produced  accurate  results  for  oil  four  points  originally 
navigated  by  the  SSEO  and  accurate  line  number  solutions  for  TA.MU 
coordinates  in  a  later  chock.  A  check  of  the  element  coordinate  for 
TAMU  was  not  accurate  in  the  later  chock  duo  to  loss  of  the  original 
input  "gamma  shifts"  used  to  produce  the  SSFC  solutions. 

The  procedure  adopted  in  the  remainder  of  this  study  depended 
upon  the  following  assumptions: 

1)  The  navigation  solutions  produced  for  the  four  points  by 
the  SSFC  are  accurate .  Young  (1970)  suggests  an  accuracy  to  within 
one  lino  number  and  within  two  element  numbers  in  tire  full-resolution 
visible  data.  The  numbering  system  for  lines  and  elements  in  the 
infrared  data  is  based  on  that  for  the  visible  data  and  is  incremented 
by  a  factor  of  four  which  suggests  an  accuracy  to  the  nearest  pixel 
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in  t ho  inn urod  data, 

2)  Intel potation  in  time  using  a  fourth-order  polynomial  pro¬ 
ducer  accurate  .solutions  for  intermediate  image  times.  This  assump¬ 
tion  is  required  as  a  compromise  to  avoid  the  computer  costs 
associated  with  providing  8i>'£C  solutions  for  images  at  15-min 
intervals  over  A  h  and  ior  four  points  in  each  image.  The 
variation  with  time  of  the  original  SSIlG  solutions  suggests  that  this 
assumption  is  valid, 

3)  Mapping  equations  of  the  form  of  Eq .  (A)  and  Eq .  (3)  pro¬ 
vide  an  accuralv;  coordinate  t.i ansfomation  between  the  line  end 
element  coordinates  of  the  satellite  image  ai ray  and  the  (>:,  y) 
coordinates  of  the  digital  radar  grid.  MAGMA  P  3  mapping  seditious 
provided  accurate;  results  for  the  SSEC  solutions  for  all  lour  navigated 
points  and  for  TAMU  coordinator .  Comparisons  of  contoured  fields 

of  infrared  data  mapped  using  MAC  MAP  3  solutions  to  hard  coj  •'/ 
enhanced  GOES  infiorod  imagery  also  confirm  tiro  general  accuracy 
of  the  mapping  procedure.  It  is  evident,  however.,  that  the  earth 
curvature  considerations  used  in  establishing  the  radar  grid  for 
MAC  MAP  3  would  not  compensate  for  such  effects  at  distances  far 
removed  from  the  center  of  the  grid.  Tor  studies  of  data  within  200 
km  of  TAMU  the  MACMAP  3  procedure  appears  to  yield  satisfactory 
results. 


Actual  production  of  maps  of  satellite  data  depended  upon 
development  of  programs  designed  to  use  MAC  MAP  3  solutions  for 
the  coefficients  of  Eq.  (4)  and  Eq.  (5).  The  original  program  was 
called  MACVIS  and  was  designed  to  map  full-resolution  visible  GOES 
data  to  the  radar  grid.  Problems  with  brightness  normalization  due 
to  a  very  large  solar  zenith  angle  led  to  abandoning  the  use  of  th.: 
visible  data  after  MACVIS  was  developed.  The  MACIR  computer  pro¬ 
gram  was  dc '/'..‘loped  from  the  MACVIS  program  and  is  designed  to-use 
solutions  to  Eq.  (4)  and  Eq.  (5)  to  map  GOES  infrared  satellite  data 
into  the  radar  grid  centered  on  TAMU. 

MACIR  is  a  modular  computer  program  using  two  function  sub  ¬ 
programs  to  produce  solutions  to  Eq.  (4)  and  Eq.  (S)  and  using 
subroutines  to  save  an  array  of  satellite  data  that  just  overlaps  the 
area  of  the  radar  grid  that  is  of  interest  and  to  map  the  saved  data 
array  into  the  radar  grid.  Subroutine  POUND  c  atormines  bounding 
line  and  element  numbers  so  that  the  saved  satellite  data  can  bo 
extracted  from  the  image  sector  records  on  the  IviclDAS  image  save 
tapes.  Subroutine  DATSAV  reads  the  image  sector  records  and  saves 
the  GOES  data  between  the  bounding  line  and  element  numbers 
determined  by  BOUND.  Subroutine  DATMAP  computes  a  line  and 
element  number  for  each  (x,  y)  point  in  the  portion  of  the  radar  grid 
selected  for  analysis  and  assigns  the  value  of  the  pixel  corresponding 
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to  that  (line,  clement)  pah  to  tlic  (>:,  y)  point  ir.  the  radar  gHd. 

Tho  data  thus  mapped  arc  stored  in  on  array  called  IGOI'S .  Ssbrou- 
tine  OUT PUi  is  used  to  print  maps  of  the  TOOL'S  array  and  is  readily 
adapted  to  varying  grid  spacing. 


The  MAC1R  program  was  used  to  produce  maps  of  the  003  3  GMT, 
0030  GMT,  0313  CMT,  and  0230  GMT  COES-  hast  images  for  3  May 


1978  for  this  investigation .  Maps  of  all  the  images  between  0015 


GMT  and  03-13  GMT  also  were  produced  lor  a  diff- 
maps  consisted  of  point  values  of  GOES -East  infr 
2-km  intervals  for  this  investigation  and  4-km  inv- 
maps.  A  listing  of  the  MAC  MAP  3  cud  MACIR  pre 


rent  study.  The 
ned  data  printed  at 
reals  in  the  ether 
ji'i'SiS  is  contained 


in  Appendix  A. 

c.  Correction  for  Displacement  Error  Duo  to  Cloud  Height 

Cloud  interpretation  near  tho  horizon  of  a  geostationary  satellite 
imago  is  comi>licatcd  by  the  effect  of  curvature  of  the  earth  surface. 
Near  the  satellite  subpoint  on  the  earth's  surface,  the  apparent  and 
actual  cloud  positions  coincide.  However,  the  difference  between 
the  apparent  and  actual  positions  increases  as  the  cloud's  distance 
from  the  subpoint  increases.  In  order  1o  locate  accurately  a  cloud 
feature  in  a  geostationary  satellite  image  it  is  necessary  to  com¬ 


pensate  for  this  effect  by  making  a  displacement  correction  toward 
the  satellite  subpoint.  Weiss  (1978)  gives  a  method  of  correcting 


SMS  images  foi  this  effect  and  that  method  was  used  to  compute 
displacement  errors  lor  various  cJoud  heights  using  a  GOES -Sant 
subpoint  of  0.335°  S,  75.021°  \V .  Displacement  corrections  were 
computed  for  TAMU  for  various  cloud  heights  and  arc  listed  in  Table 

1. 

Table  1 ,  Di  spin  cement  cor recti  on:;  for  cloud  h?  jc  h1_ , _ 


Cloud  Height  (km) 

Displacement  Collection  (k.rn) 

10 

9.17 

12 

11.00 

14 

12.83 

16 

14.66 

Cloud  height  was  determined  by  converting  an  average  digital 
count  for  each  image  to  an  equivalent  black  body  temperature  and 
selecting  the  ICAO  stondaid  height  corresponding  to  that  temp  Taturc 
in  the  2350  GMT  College  Station  sounding.  1'or  temperatures  near 
and  above  the  tropopause  a  different  pioccdurc  wa s  used.  The 
saturation  adiabat  from  the  convective  condensation  level  was  used 
to  determine  an  equilibrium  level  on  the  sounding.  Above  this  level 
a  subjective  estimation  of  the  temperature  of  convective  tops  bas^d 
on  approximately  20  per  cent  mixing  with  environmental  air  was 


used  to  determine  a  temperature -height  relationship. 


Rather  than  attempt  to  shift  each  point  in  each  satellite  imago 
toward  the  satellite  subpoint  an  amount  dependent  upon  the  indivi¬ 
dual  digital  count  at  each  point,  a  uniform  shift  of  the  elate,  toward 
the  satellite  subpoint  based  on  the  average  digital  count  for  the 
image  was  used.  This  procedure  was  adopted  to  avoid  computations 
problems  and  to  prevent  creation  of  "shadows"  in  highly  contoured 
infrared  data  fields.  Table  2  lists  the  magnitudes  of  the  corrections 
applied  to  the  images  in  this  investigation. 

_ Table  2.  Displacement  correction  for  each  ivn-uy. _ 


Image  Time  Average  Digital  Count  x  y 

Displacement  Di  spin  com  art 
(hm)  (km) 


001 5Z 

190 

5 

8 

0030Z 

204 

G 

9 

0215Z 

211 

6 

9 

02  3  0Z 

218 

7 

11 

After  the  GOES  infrared  images  were  mapped  to  the  radar  grid 
and  corrected  for  displacement  error  due  to  cloud  height,  the  data 
for  these  images  were  tabulated  along  with  corresponding  radar  data 
in  the  three  radar  data  sets.  These  data  then  formed  the  three 


data  sets  examined  in  this  investigation. 
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CHAPTER  IV 

COMPARISON  OF  DIGITAL  RADAR  ADD  SATELLITE  DATA 

Comparison  of  the  field?;  of  digital  radar  and  digital  satellite 
data  was  performed  by  first  qualitatively  comparing  contoured  maps 
of  the  data  and  then  objectively  comparing  the  fields  by  computing 
correlation  coefficients  between  the  radar  and  satellite  data  using 
tho  SAS  76  statistical  analysis  system  (Barr  cl  al . ,  1  97  0)  as  imple¬ 
mented  on  the  Amdahl  470  V/(>  computer  at  TAMl). 

1.  The  Pa  sis  for  the  Comparison 

Meteorological  radars  examine  the  structure  of  weather  systems 
by  emitting  electromagnetic  energy  in  the  microwave  portion  of  the 
electromagnetic  spectrum  and  measuring  the  amount  of  that  energy 
scattered  back  to  tire  radar  by  hydrometeors  of  various  sizes  and 
shapes  associated  with  the  weather  system.  Meteorological 
satellites  examine  weather  systems  by  sensing  the  intensity  of  scat¬ 
tered  or  emitted  electromagnetic  energy  from  the  cloud  structures 
and,  ultimately,  from  the  hydrometeors  that  make  up  the  clouds 
themselves.  Sensors  aboard  meteorological  satellites  typically 
respond  to  electromagnetic  radiation  in  the  visible  and  infrared 
portions  of  the  electromagnetic  spectrum  although  ultraviolet  and 
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microwave  son  son;  have  boon  used  for  special  applications. 

Since  both  meteorological  radars  and  meteorological  satellites 
examine  weather  systems  by  the  use  of  scattered  or  emitted  electro¬ 
magnetic  radiation  at  different  wavelengths, it  is  reasonable  to  ex¬ 
pect  that  some  relationship  should  exist  between  meteorological 
radar  data  and  meteorological  satellite  data  that  are  coincident  in 
time  and  space.  The  basic  purpose  of  this  investigation  is  to  deter- 

A 

mine  if  such  a  relationship  exists  and,  if  so,  to  measure  the  strength 
of  such  a  relationship. 

To  examine  the  relationship  between  meteorological  radar  data 
and  meteorological  satellite  data,  two  tilt -sequences  of  digital 
radar  data  and  four  infrared  GOES- Hast  satellite  .images  were  exa¬ 
mined.  The  first  tilt-sequence  was  selected  early  in  the  data 
collection  period  when  radar-detected  weather  activity  in  the  area 
was  relatively  weak  and  just  beginning  to  develop.  GOES -II a st 
images  collected  just  prior  to  initiating  the  tilt-sequence  and  at  its 
conclusion  wore  used  for  the  comparison.  The  second  tilt-sequence 
was  selected  from  a  later  time  interval  in  which  line  convective 
activity  in  the  TAMU  area  had  reached  near-maximum  intensity  for 
the  data  collection  period.  Comparison  GOES -East  infrared  images 
also  were  selected  for  times  closely  corresponding  to  initiation  and 
termination  of  the  tilt-sequence.  The  results  of  the  comparisons  of 
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these  tilt-sequc.nces  with  the  corresponding  satellite  data  are  pre¬ 
sented  in  the  following  sections. 

2.  First-Tilt- sequence  and  Corresponding  GOES  Images 

The  first  digital  radar  tilt-sequence  was  collected  from  0022 
GMT  to  0030  GMT  on  3  May  1978.  The  data  wore  collected  in  14 
1°  inclinations  of  the  radar  antenna  with  each  change  in  antenna 
angle  being  made  at  the  conclusion  of  a  360°  swoop.  Data  from 
the  northwest  sector  of  the  tilt-sequence  were  processed  by  the 
Sieland  (1977)  computer  program.  The  0015  GMT  and  0030  GMT 

I 

GOES-East  infrared  images  wore  used  fox  comparison  with  this  tilt— 
sequence.  Maps  of  zero-tilt  reflect jvily  (dBZ) ,  low-,  mid-,  and  • 
upper-level  FVSZ  (unZ) ,  and  VI L  (kg  m  )  are  presented  in  Figs. 

3,  4,  5,  6,  and  7,  respectively.  Contoured  maps  of  infrared 
digital  count  for  the  0015  GMT  and  0030  GMT  GOES-East  images 

for  corresponding  areas  are  presented  in  Figs.  8  and  9,  respectively. 

Conversion  from  infrared  digital  count  to  equivalent  black-body  tem¬ 
perature  can  be  approximated  by  using  Fig.  10. 

Significant  features  to  note  in  the  digital  radar  data  include  a 
marked  absence  of  indications  of  echo  activity  in  the  upper-level 
PVSZ  map  (Fig.  6)  and  the  VIL  map  (Fig.  7).  Echo  activity  of 
weak  to  moderate  intensity  is  indicated  in  the  zero-tilt  reflectivity 
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I’ig.  G.  Tilt-sequence  one,  upper-level  PVSZ  (dBZ).  Date  collec 
on  3  May  1978  from  0022  GMT  to  0030  GMT. 
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map  (lie.  3),  the  low-  level  PVS/1  ::j;.p  (l‘jy .  -1),  nnd  in  tfv  ::i id -lc v«  1 
l’VSZ  map  (Fig .  5).  The  bost-o.'g :.ni  ::.od  coho  activity  is  l."oa!:cd  In 
the  western  and  soulh-coutrul  portions  of  1' i<_> .  3  anti  -l. 

The  001  3Z  infrared  image  (Fin.  0)  ha  a  a  n.t  ar-uniform  or  .:dient 
of  digital  count  from  high  value.;  to  low  from  r.enhwcst  m  southeast. 
The  sharp  int  rouse  in  gradient  in  the  couth;-:: stem  portion  of  Lise 
imago  is  probably  caused  by  the  censor  "loch. in;:  at"  lower- level 
cloud  or  the  side  of  convective  at.  1  ivity .  The  shapes  of  the  contours 
for  digital  counts  of  2 OF.  and  310  r-vgg act  !.am:  acsoci ..tion  with  con- 
voctive  activity  as  evidenced  in  that  area  by  Fly  a.  3,  4,  and  5. 

The  003013  iiifran  image  {!  jg  .  3)  lias  u  me  live  jy  ih:t  yruwiont 
of  digital  count  over  much  of  the  snapping .  The  closed  210  contour 
suggests  its  relation  to  convective  cells  shown  in  the  same  area  in 
Figs.  3,  A,  and  5.  The  relative  flatness  of  the  gradient  of  digital 
counts  on  the  order  of  200  and  the  disappearance  of  the  low  values  of 
digital  couni  in  the  southeast  corner  of  the  image  suggest  the  spread¬ 
ing  of  cirrus  cloud  cover  over  much  of  the  area.  Obscuration  by  the 
cirrus  prevents  direct  observation  of  low-  and  mid-level  convective 
activity  by  the  infrared  sensor.  This  effect,  when  coupled  with  the 
coarser  resolution  of  the  infrared  sensor  as  compared  with  the  digital 
radar,  is  considered  to  be  one  of  the  chief  problems  in  attempting  to 


relate  the  two  types  of  data. 
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3.  Tenth  Tilt- sequence  and  Corrc spon  ding  GOES  lieges 

The  tent])  digital  radar  tilt- sequence  was  collected  from  0216 
GMT  tc»  0233  GMT  on  3  May  1978  using  14  1°  changes  in 
antenna  inclination.  Data  from  the  northwest  and  southwest  sectors 
of  the  tilt-sequent;-:;  were  processed  by  the  Siclend  computer  program . 
The  0215  GMT  and  0230  GMT  GOES-East  infrared  images  were  us-,  d 
for  comparison  with  both  sectors  from  this  tilt-scquer.ce.  Maps  of 
zero-tilt  reflective  /  (dRZ),  low-,  mid-,  arid  high-level  PVSZ  (dBS), 
and  V]L(!.g  m-^)  arc  pro  sen  led.  in  Figs.  13  ,  12,  13,  14,  and  15, 
respectively,  for  the  southwest  sector  a.M  in  Figs.  18,  19,  20,  21  , 
and  22,  respectively,  for  the  rorlhwer.t  sect  or .  Maps  of  contoured 
infrared  digital  count  from  the  0215  GMT  and  0230  GMT  GOES-East 
images  for  corresponding  areas  ore  presented  for  the  southwest  sector 
in  Figs.  16  and  17  and  for  the  northwest  sector  in  Figs.  23  and  ?4. 
a.  Southwest  Sector 

Significant  features  to  note  in  the  digital  radar  data  include  the 
presence  of  strong  line  radar  echoes  to  the  south- sout invest  of  T/TMU 
in  all  of  the  digital  radar  maps.  The  two  areas  of  heaviest  convective 
activity  are  especially  evident  in  the  mid-  and  upper-level  PVSZ  maps 
(Figs.  13  and  14)  and  in  the  VII  map  (Fig.  15).  The  two  strong 
cells  exhibit  very  little  "tilt"  in  the  vertical  in  a  comparison  of 
Figs.  12  through  14.  Maxima  in  the  VJL  map  also  correspond  well 


Fig.  15.  Tilt-sequence  ton,  southwor;t  socto",  VIL  (kg  rn 
Data  collected  on  3  May  1978  from  021  G  GMT  to  0233  GMT. 


I’ig.  1G.  0215  GMT,  3  May  1978  infrared  image 
southwest  sector. 
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with  the  location  cf  reflectivity  maxima  in  those  fanitv  figure  a . 

High  vo lues  of  infrared  digital  court  in  both  the  021 T*  GMT  r:;-,d 
0230  Gl'.ii'  images  (Tigs.  ]  G  and  17)  lie  conic  10  hm  to  the  wost-north- 
wcst  of  radar  activity  indicated  in  Figs.  3  1  'through  15  .  The  increase  in 
maximum  digital  count  in  going  .from  the  0215  GMT  image  to  the  0230 
Glvl’i  imago  is  indicative  of  the  strong  upper -level  out  flow  associated 
with  the  line?  of  convective  activity.  Regions  of  maximum  digital 
count  in  both  images  lie  overlap  the  maxima  in  tire  radar  fields, 
especially  in  the  later  infrared  image  (Fig.  3  7). 


b.  Northwest  Sector 


Strong  echo  activity  is  also  evident  in  tin  racer  de; : 


*  f  n-~  i  lv 


northwest  sector.  Strong  echoes  are  present  to  the  north  of  17*  Mil 
with  colls  in  the  lino  activity  becoming  more  evident  in  the  mid-  and 
upper-level  i'VSZ  mops  (Figs.  2.0  and  23).  Two  cells  stand  out 
very  well  in  the  V3L  mop  (Fig.  22).  Very  lil Lie*  tilt  is  evident  in 
the  cells  to  the  north  of  T/sMl)  although  the  area  of  reflc  etivity 
associated  with  each  cell  in  Figs.  3  9  through  21  is  bros dor  with 
merging  of  the  colls  nppnicntJy  occurring  at  all  levels.  The  VIL  map 
(Fig.  22)  exhibits  two  sharply  bounded  cells  almost  due  north  of 
TAMU. 

Maximum  values  of  digital  count  in  the  0215  GMT  infrared  image 
(Fig.  23)  lie  10-15  i.mto  the  west  of  maxima  in  the  radar  maps. 


Maximum  values  of  digital  count  in  t!  Ci2?,P.  CUT  infrared  imago  Ho 
slightly  to  the  west  of  maxima  in  the  dig  Hal  radar  fields.  The 
increase  in  maximum  digital  count  from  0.'ll.r>  GUT  to  0?1C  GUT  again 
indicates  the  strong  upper -level  outflow  associated  with  the  con¬ 
vective  line. 


4.  Correlations  Between  the  Under  and  Sul  . Tide  .Oulu 

Correlation  analysis  provides  a  measure  of  the  6'  ye  n  of  asso¬ 
ciation  between  variables.  The  correlation  used  as  a  nr  '.sure  of 
association  in  the  remainder  of  this  investigation  is  IV:  iron's  pro¬ 
duct-moment  correlation  with  the  correlation  cccfioie.d  .  r.  being 


given  by 
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where  the  (>h,  Yj)  represent  paired  observations  of  the  variables,  X 
and  Y  represent  arithmetic  means  of  the  variables,  and  N  is  the 
total  number  of  paired  observations  of  X  and  Y.  The  v.duo  of  the 
coefficient,  r,  is  a  measure  of  the  degree  oT  association  of  the  vari¬ 
ables  with  a  perfect  positive  association  indicated  by  r  -  1,  a  perfect 
negative  correlation  indicatedby  r  -  -1,  and  a  complete  lack  of 


correlation  between  the  variables  indicated  by  r  =  0. 
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Correlation  coefficients  wore  computed  using  procedures  imp!':  ■ 
merited  by  the  SAS  76  statistical  analysis-  system  which  also  re!  m and 
significance  probabilities  of  the  correlation  coefficients.  The  signi¬ 
ficance  probability,  p,  of  a  correlation  coefficient  is  the  probability 
that  a  value  of  the  correlation  coefficient  nr  large  or  larger  in  abso¬ 
lute  value  than  the  one  calcttl;..*  .d  would  have  ah-,  n  by  chance,  wmo 
the  tv.e  random  v.ninMos  truly  unooiroMb -b  (berry.  _-h  ,  197C) .  Cur- 
rolatic.n  co<  fficicnls  were  ca!cul*t<  d  foi  three  bo  si-.:  data  sets:  rhe 
first  t  jU-sequenc-.  and  associaiod  GC: ..u-s,  11.:'  tenth  tilt- 
quence ,  sc. it  Invent  sector,  arid  as  social-  d  COLS  >rr-ay«  s,  and  L.h 
tenth  t i It -•  s-.. c j a •  nc.  ,  noli, in. • 'ol  ror ,  o-.  .  j.-  e-jo.1..-.-.'...  rmayrs.  / - j . 
values  of  r  >,voie  comput  -1  ushsi  point  v;  !  of  the  rs.de  r  and 
satellite  variable  s . 

An  adddionnl  sebdliio  rial  •  v.u'iabh  was  artifie.l’y  cr.  ati  d  h 
each  data  set  by  subtracting  i:<  hi;  value:'  of  tlio  digit--.!  '•'oust  of  b 
later  GO:', 6  mage  from  con  os.  on-ling  volu-as  cf  dieiu-l  conn*  in  i  i 
earlier  COLS  image.  Correlation  coefficient  a  also  v.vr-  -  computed 
between  this  differ. -nee  in  digital  count,  d,  and  the  digital  radar  data . 
a.  First  Tilt-sequence  and  Corresponding  GOl'f!  Images 

The  results  of  the  correlations  of  the  0015  GMT  infrared  data, 
the  0030  GMT  infrared  data,  and  d  with  corn,  spending  digital  radar 
data  arc  tabulated  in  Table  3.  The  best  correlations  wore  found 
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between  the  0015  C-MT  end  0030  GMT  image?,  and  mid- level  PVSZ  data 
Weak  and  sometimes  negative  correlation  coefficients  were  found  be¬ 
tween  both  infrared  imago;,  and  zero --tilt  reflectivity  and  low -level 
PVSZ.  V/euk  positive  coi  re  lotions  are  indicated  for  both  images  and 
upper-  leve  l  PVSZ  and  VII,.  The  results  for  the  point  differences  in 
digital  counts  between  images  and  the  rad, a  date  arc  ambiguous  with 
a  weak  positive  correlation  between  d  and  zero-tilt  reflectivity  end 
low-level  PVSZ  and  negative  correlations  with  t!v  remaining  radar 
data  fields.. 

b.  Tenth  Till- sequence  ,  Southwest  Sector,  and  Cora  spanning  GOTO 
Image s 


The  rer 

ulls  c.f  cora. lei  ions  of  tie  hi 

'•  at ; 

’  i  nfra’ .  k  .1 

,  tin 

0230  GMT  i 

idrarod  data,  and  d  with  the 

oa'.-:g 

••■•diag  digit.. 

1  X.du 

data  from  LI 

!  southwest  sector  of  th<  te’ 

,T.  till 

:V  jiienc  •  ex 

i-be 

in  Table  -1.  The  best  corn  lotions  in  this  dot  j  t  v, ,  rc-  !.  etw  ■  i  mid 
level  and  upper-level  1’VbZ  and  u,  the  JLffex  nci;  between  digital 
counts  for  tie,  0330  GMT  and  021  5  CM. T  image;;  .  Voir  porilivi  cor¬ 
relations  also  were;  found  b  tween  r!  and  VI 1-  and  between  the  0230 
GMT  image  and  zero-tilt  reflectivity  and  low-level  PVSZ.  Somewhat 
weaker  correlations  we  re  found  between  the  0330  GMT  image  and  mid- 
and  upper-level  PVSZ  and  VIL.  Weak  positive  correlations  woic  found 
between  the  0215  GMT  image  and  zoio-tilt  reflectivity  and  low-level 


To  bio  4.  Values  of  r  and  p  for  tl-c*  tenth  tilv-s-'quc-neo 
and  corrof.pc-ndiiHj  COLS  irsidyos  (southwest 
_ recto; ) . _ _ _ _ _ _ _ _ 
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Radar  Data 
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0 . 0CO1 

Upper- -love!  PVSZ 

0.G2-: 

0. 00,01 

VJL 

0.  5  2  7 

0.(  oOj 

Numb'-r  of  GOrer  ml  ion:.;  ■  1900 


PVSZ.  Negative  con di Lions  misled  between  the  0215  GMT  iir.o:.;-..- 
mid-level  and  upper- level  FVSZ  iind  Vjf,.  Zero -till  n  \f  iectivity  on-! 
low- level  PVSZ  were  essentially  uncorrelnted  with  d. 
c.  Tenl.ii  Tilt- sequence,  Northwest  Secioi  ,  and  Corresponding  (20 
Images 

Table  5  lists  the  result:.  or  eorre  lei  ion:-  of  023  5  C'.-iV  mfree  2 
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proved 
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in  examining  the  structure  of  convective  stor 

;n  systems 

,  and  iii  >: 

mating  the  amount  of  precipitation  from  convective  storm  system: 

'la bio  5.  Values  of  r  and  p  lor  tbo  tenth  tiH-rvquence 
and  corresponding  COPS  images  Ciour.  west 
_ sector) .  _ 


Satellite  Data 

Radar  Data 
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0215  C lvIT  Image 

Zero  -tjit  refit  ctivity 
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Low-!’  \n  ]  RVSZ 
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Eniv  nrirg  i  1  > f j  :•  >  _  cl  Images  requires  b.-e ting  v;  rj<  -us  "thr>  .•  3 :■; > !tl:: "  ei 
the  data  and  assigning  vuiious  shades  of  grey  in  the  "9 re y  sen)-.  - "  to 
data  bounded  !.>y  tho  thresholds.  These?  thresholds  may  the.:  bo  iu;<  d 
to  pi  educe  patterns  in  tho  infrared  data  cherocLci  istio  of  vatious 
stages  in  the  d'  .velopmeni  of  jrhonomen.'i  such  ng  hna  icanc;:  and  squu  •  ’ 
liif.'  j  . 


An  example,  of  an  enhancement  aj.pl  lets  to  in  ft  overt  imwery  ir,  the 
MB  enhancement  curve  n.ho.vn  in  Fin.  2b  and  tubulated  in  Tuia-  „■  6, 
The  MB  onh  irxvment  is  i.ppUed  to  COiib  infrared  imagery  or  a  men;; 
of  delineating  areas  of  thunderstorm  activity  in  the  images.  The  wry 
cold  tens  tvoiwliv  tir.sooic.led  with  «n  eas  oi  intense;  cotivecfu  "  net  i- 


vity  arc  ti.r.-r Folded  in  f !.  •  grey  scale  areas  of  dai 5;  and  eh  1. 
seen  against  the  lighi<  r  har). ground  d  lower-level, warmer  cicada 
the  hard  copy  imagery.  A  re  vi  real  of  t li- •  rcale  from  black  to  whit, 
used  to  identify  the  Stroup  !y  pencil  .dive  lops  indie. itivo  of  sev<  re 
thunderstorm  activity  (Ccrbail  ft  a] .  ,  13/0). 
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a.  Thresholding  t he  infinite.  Data 

Selected  digital  count  thresholds  correr- ponding  tc  levels  em¬ 
ployed  in  tho  MB  enhancement  curve  wore  ajiplied  to  the  data  used 
in  this  invest iqation .  All  infrared  data  having  digital  counts  less 


than  the  threshold  digital  count  and  corresponding  radar  dal  a  were 


exclude  d  from  that  step  in  the  c-unlysir 
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other  factors.  Thresholding  the  0030  GMT  result"  produced  siin 
results  for  the  193  digital  cour.l  Ihi oshold  ir>  that  all  cm  re1  a  Lion 
were  somewhat  improved.  The  201  and  20G  digital  count  thror.ht 
improved  correlations  with  zero-tilt  reflectivity  and  low-'1'  vel 
while  the  correlation  with  mid-level  PVCZ  diminished.  Ti.  ran 
with  upper-level  PVGZ  remained  orrenthdly 
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the  thresholding  process.  The  conflation  wi ‘ !  VM.  was  generally 
improved  by  Uio  thresholding  proem.;;  with  tb.-  heel  coue’niion  -,t 
digital  count;;  at  and  above  298. 
b.  Thresholding  the  Digital  Ruder  Data 

The  tluoshclding  piccess  was  extended  to  the  digital  uidar  data 
as  well  and  con  elation  co--rTici  out;,  were  co  -puled  b’-tw:  a  the  radar 


and  satellite  data.  The  initial  thresholding  v 
mining  only  radar  data  with  nonzero  value-.*;  fan 
few  isolated  coses  a  slight  i:npro\  esv.  t  in  co 
satellite  dot.)  and  zero-tilt  ’.->flcot,'vity  was  o? 
cases  l.ne  coiiolation  was  diminish..  .1,  mark.:;. 


’■  -  p.-ronn.-.-d  by  exu- 
the  Vurir.bk's.  In  a 
lotions  between  the 
served .  In  all  other 


do  In  ro’.'ri  f n rl ivncb’ c 


To  examine  further  tin  effect  of 
data,  arbitrary  threshold;;  woo  setae 
computed  between  the  radar  and  net: 


hr-.shohnng  Die  digital  radar 
.;:d  and  c.  v  r.  la  lion  coefficients 
iiledato.  Threshold  levels 


of  1  0  dBZ  were  selected  for  zero-tilt  reflectivity.  Jew -level  PVSZ, 


and  mid- level  PVSZ  data,  h  thresliold  level  of  5  c’liZ  was  selected 
for  the  upper  level  PVSZ  data  ana  of  5  kg  m  for  the  VIL  data  . 
Correlation  coefficients  computed  for  the  threshoklod  radar  data  and 
corresponding  infrared  images  are  listed  in  Table  8. 

Thresholding  the  radar  data  in  the  first  tilt-soqucnce  improved 
the  correlation  between  zero-tilt  reflectivity ,  low-Iovel  PVSZ,  and 
mid-level  PVSZ  and  both  infrared  images.  No  improvement  was 


t  \ 
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obt>.:i nod  in  the  other  correlations  or  feu  correlation  r.  of  the  ruder  data 


with  d. 


Thresholding  the  radar  data  in  tin  southwest  sector  o.f  the  tenth 
tilt-  sequence  improved  the  conv  latkw  cf  d  with  y.<  reflectivity 

and  low-loved  PVSZ  but  did  not  improve  any  other  correlation  in  the 
data  sot.  Thresholding  the  radar  d..!a  in  the  northwest  seder  of  the 
tilt-sc  (iuctu-e  .irurov'.d  the  cosrelatien  o£  both  infra:  ed  in*,  eg  os  with 
zero-tilt  refh •c.Tiviiy  and  of  the  0218  GMT  image  with  low-level 
PVSZ.  The  remaining  co;  relations  w._re  not  imp.e.  wed  by  the 


thresholding  prr  ocss. 

Simultaneous  thiv  nhokiiny  of  beta  the  rijnita.l  r.idar  and  bigUa! 
satellite  data  was  not  oth  moled  sit  re  only  data  from  the  first  tilt — 
sequence  and  corresponding  imagery  could  have  born  used.  Reduced 
sample  size  in  such  thresholding  p>-c.;  ably  would  Sieve  resulted  in 
misleading  or  erroneous  results,  particularly  for  comparisons  willi 
upper-level  PVSZ  and  VIL. 


Cl  3  A  ITER  V 


CONCLUSION  3  AMD  RE  CO  Tv; 


DAT  JONS 


1.  Com:!  visions 

The  primary  objective  of  tins  Investigation  was  to  compare  digital 
meteorological  radar  dotn  will-*  infrared  imagery  Xro:r>  a  men.  oroiogicn  l 
satellite  in  geostationary  orbit  by  determining  tlic  con  o’.dion 
between  point  values  of  infrared  digital  count  and  point  values  of 
radar  reflectivity  and  vertically  .integrated  liquid  wafer  content  (VIE). 
Digital  do  fa  from  the  10  cm  radar  of  the  Dupwimoni  of  LVioorclogj  at 
T/'iIviTJ  were  processed  V>v  ri  niorliTir  f]  vorr. ■io:*?  .if  (1  ^ 77)  < 

put  or  program  and  displayed  as  output  maps  of  point  values  cl  zero- 
tilt  reflectivity,  low-,  mid-,  and  upper-level  rVf’Z,  and  VI L,  Digital 
GOES- East  infrared  satellite  data  from  I.fcIAuS  image  save.-  tapes 
processed  to  produce  maps  of  point  values  of  infrared  digital,  count 
mapped  to  the  same  projection  and  giid  interred  used  for  the  digital 
radar  data.  The  maps  of  digital  radar  and  digital  satellite  images 
formed  the  basis  for  comparisons  between  the  two  types  of  data. 

This  investigation  has  led  to  the  following  conclusions: 

1)  Large-scale  features  evident  in  the  digital  radar  data  arc 
also  apparent  in  the  infrared  imagery.  The  obscuring  of  low  and 
middle  level  convection  by  cirrus  overcast  and  the  different 


ro«t>J lit in  the*  imrurod  GOES  imagery  (—-8  km)  end  the  di  glial 
radar  (— 1  km  interpolated  to  a  2-kin  grid  spacing)  make  areas  of  wo 
convection  that  arc  apparent  in  the  radar  data  difficult  or  impossible 
to  locate  in  the  infrared  imagery.  This  p: oblom  is  probably  me- a*  1  rr- 
portent  to  rainfall-ost imutiou  efforts  Hum  lo  efforts  to  locate  arid  •;.i 
sure  the  intensity  of  seven  convocl ivo  activity.  Tie-  fair  corral;. tie 
between  mid-level  PVSZ  and  the  infrared  do  to  in  such  am  instance 
may  provide  a  v.my  to  partially  eliminate  ties  problem,  particularly 
when  the  infrared  data  ore  enhanced  by  thi'.mboklir.e . 

2)  Although  only  one  set  of  2001  observations  could  be  arm 
mined,  thresholding  the  inhaied  data  significantly  improved  the 
correlation  botwc  on  the  ladar  and  saloliic.  d  .in ,  especially  the; 
correhii: ion  with  :a  lo-tilt  reflectivity  and  h/.v-level  PV SZ .  Chong 
Rod  or:  huh:  (1977)  compared  NO.hr.  2  infra! cd  imagery  wit],  radar  rain¬ 
fall  ami  determined  that  thresholding  the  satellite  imagery  did  r.ol 
improve*  ‘ho  correl  iLion.  The  results  of  this  investigation  suggest 
that  thresholding  the  infrared  GOES  imagery  should  improve  the  (cor¬ 
relation  between  the  satellite  data  and  radar  rainfall.  The  different' 


in  results  may  be  due  to  procedural  differences  or  to  differences  in 
the  satellite  sensors  since  the  infrared  sensor  aboard  the  NOAA  2 


satellite  has  a  higher  resolution  than  the  GOES  VISSR. 

3)  Com  lations  between  tho  radar  and  satellite  data,  while 
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Rod  Henderson  is  married  to  the  former  Sandro  hoggins  of 


